arXiv:1504.04028v2 [astro-ph.SR] 30 May 2015 


Time variation of Kepler transits induced by stellar spots — 
a way to distinguish between prograde and retrograde motion. 

II. Application to KOIs 

Tomer Holczer^, Avi Shporer^’^’^°, Tsevi Mazeh^’^, Dan Fabrycky®, Gil Nachmani^, Amy 
McQuillan^, Roberto Sanchis-Ojeda®’^°, Jerome A. Orosz^, William F. Welsli^, Eric 

B. Ford®’®, and Daniel Jontof-Hutter® 

Received_; accepted_ 


^ School of Physics and Astronomy, Raymond and Beverly Sackler Faculty of Exact 
Sciences, Tel Aviv University, Tel Aviv 69978, Israel 

^Division of Geological and Planetary Sciences, Galifornia Institute of Technology, 
Pasadena, GA 91125, USA 

®Jet Propulsion Laboratory, Galifornia Institute of Technology, 4800 Oak Grove Drive, 
Pasadena, GA 91109, USA 

'^Kavli Institute for Theoretical Physics, University of Galifornia, Santa Barbara 
^Department of Astronomy and Astrophysics, University of Ghicago, 5640 South Ellis 
Avenue, Ghicago, IL 60637, USA 

^Department of Astronomy, University of Galifornia Berkeley, Berkeley, GA 94720, USA 
^Astronomy Department, San Diego State University, San Diego, GA 92182, USA 
®Department of Astronomy and Astrophysics, The Pennsylvania State University, 525 
Davey Laboratory, University Park, PA 16802, USA 

®Genter for Exoplanets and Habitable Worlds, The Pennsylvania State University, Uni¬ 
versity Park, PA 16802, USA 
Sagan Fellow 





2 


ABSTRACT 

Mazeh, Holczer, and Shporer (2015) have presented an approach that can, 
in principle, use the derived transit timing variation (TTV) of some transiting 
planets observed by the Kepler mission to distinguish between prograde and 
retrograde motion of their orbits with respect to their parent stars’ rotation. 
The approach utilizes TTVs induced by spot-crossing events that occur when the 
planet moves across a spot on the stellar surface, looking for a correlation between 
the derived TTVs and the stellar brightness derivatives at the corresponding 
transits. This can work even in data that cannot temporally resolve the spot¬ 
crossing events themselves. Here we apply this approach to the Kepler KOIs, 
identifying nine systems where the photometric spot modulation is large enough 
and the transit timing accurate enough to allow detection of a TTV-brightness- 
derivatives correlation. Of those systems hve show highly signihcant prograde 
motion (Kepler-17b, Kepler-71b, KOI-883.01, KOI-895.01, and KOI-1074.01), 
while no system displays retrograde motion, consistent with the suggestion that 
planets orbiting cool stars have prograde motion. All hve systems have impact 
parameter 0.2 < 5 < 0.5, and all systems within that impact parameter range 
show signihcant correlation, except HAT-P-llb where the lack of a correlation 
follows its large stellar obliquity. Our search suhers from an observational bias 
against detection of high impact parameter cases, and the detected sample is 
extremely small. Nevertheless, our hndings may suggest that stellar spots, or at 
least the larger ones, tend to be located at a low stellar latitude, but not along 
the stellar equator, similar to the Sun. 

Subject headings: planetary systems — techniques: photometric — stars: spots — 


stars: rotation 
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1. Introduction 

Observed characteristics of exoplanet systems may contain clnes to their formation 
process and orbital evolntion history. One of the promising clnes is stellar obliqnity, the 
angle between the stellar spin and planetary orbital angnlar momentnm, also called the 


to resnlt in well a 


Winn et ah 


201Qa; 


igned systems (e.g., 


Dawson 


Fabrvckv fc Tremaine 


20071: 


Cresswell et al. 

2007; 

Nagasawa et al. 

2008 


2014t). and others expected to p rodnce misaligned orbits (e.g., 


Naoz et al. 


2011 


Batygin 


2 OI 2 I) . Those theories are now being 


pnt to the test ns i ng the increasing sample o f host stars with measnred obliqnity (e.g., 


Triaud et al. 

2010; 

Albrecht et al. 

2012, 

2013) 


The growing 


s ample has led to the detection of several possible observational trends. 


Fabrvckv fc WinnI (l2009l) have identihed early on that the sample of obliqnity measnrements 


is better explained as being composed of two different nnderlying distribution^ one we ll 


aligned and the second isotropic, than a single distribution. Later, 


Winn et al. 


fl2nina 


suggested that cool stars, below Tes ~ 6,250 K, reside in well aligned systems and that the 


obliquity of hotter stars cover a wider range. This was su pported by a s 
sample by 


Albrecht et al. 

(2012 

) and a statistical study by 

Mazeh et al. 

(2015b1 


trend, identihed by 


udy o f a larger 


Hebrard et al. 


fj201l[) . shows that planets in a retrograde orbit are less 


massive than a mass cut-off of about 3.5 Mj, while planets above that threshold tend to be 
in prograde orbits. This suggests that different orbital evolution p roces ses are at play above 


and below that mass threshold. More recently. 


Morton fc WinnI (12014) showed that stars 


hosting a multi planet system (as detected by Kepler) tend to be more well aligned than 
stars hosting a single planet system. 


The need for a large sample of systems with a measured host star obliquity has 
led to the development of several methods for measuring it. Those include utilizing 






































































4 


the Rossiter-McLaugh li n (RM) effect (e.g. 


Gizon fc Solank 


2003 


stellar rotation (e.g ., 


Chaplin et al 


Schlaufman 


2013 


Mazeh et ah 


20101 : 


Gaudi fc Winn! 12007 h asteroseismology 


Van Evlen et a 


Hirano et alJ 


2015b|), stellar gravity darkening flBarnes 


2012 


20 


2014 


2009 


Lund et al 


20 


Morton fc Winn 


Szabo et ah 


stellar activity in the form o f spots (e.g., 


Sanchis-Oieda fc Winn 


2012 


20131 


2011 


Dfeert et al 


2011 


e.g., 


iiu 


2014- 


Barnes et ah 


Nutzman et al. 

2011: 

Sanchis-Oieda et al. 

2011 

.2011: 

Deming et al. 

2011; 

Sanchis-Oieda et al. 


The methods above are either statistical in nature, so no conclusions can be made about 
the configuration of a specific system, or require additional data—spectroscopic and/or 
Kepler short cadence light curves. Mazeh, Holczer, and Shporer 2015 (Hereafter Paper I) 
presented a new method to distinguish between prograde and retrograde motion of transiting 
systems for photometrically active host stars. The method utilizes spot-crossing events, 
where the transiting object, a planet in our case, moves acr oss a dark spot on the sur f ace o 


the host star. Unlike the approach of 


Nutzman et al. 


see also 


Sanchis-Oieda et al. 


(2011) and 


Sanchis-Ojeda et al. 


(12011 


201211 ■ which usually requires Kepler short cadence data, this 


method can use Kepler long cadence data, available for all objects observed by Kepler. 

Briefly, the method presented in Paper I is based on measuring two parameters for each 
transit event. The first is the local stellar brightness temporal derivative at transit time, or 
simply the local slope, a parameter that captures the host star’s rotation phase as the spot 
is rotating towards or away from the center of the stellar disc. The second parameter is 
the transit timing v ariation (TTV), assumed to be induced by the spot-crossing event (e.g.. 


Oshagh et al. 


201311 . Spot crossings distort the transit light curve shape, and although with 


long cadence data they usually cannot be fully resolved, they result in a small mid-transit 
time shift when fitted with a simple transit light curve model that does not account for the 
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spot-crossing event. The TTV reflects the phase within the transit where the spot-crossing 
occurred. As shown in detail in Paper I, a negative (positive) correlation between the 
TTV and the local slope indicates a prograde (retrograde) motion, under some simplistic 
assumptions. 

Here we present the application of the method presented in Paper I to Kepler KOIs. 
In Section 2 we describe our analysis and the way we identify objects with statistically 
significant correlation between the TTVs and local slopes. In Section 3 we discuss the transit 
impact parameter of the systems expected to be sensitive to TTV due to spot-crossing. 

In Section 4 we present the six systems for which we detect a statistically significant 
correlation, while in Section 5 we discuss the systems that are expected to show significant 
TTV due to spot-crossing although no significant correlation is detected. Section 6 presents 
a brief discussion and a summary. 


2. Identification of KOIs with significant Correlation 

As mentioned above, our method uses two parameters measured for each transit 
event—the TTV and the local slope of the stellar brightness at the time of the transit. 
Section [2T] briefly describes the catalog of transit timing measurements we use here, and 
Section 12.21 describes our measurement of the local photometric slopes. In Section 12.31 we 
derive the correlation between the two sets of parameters for each KOI, and point to the 
significant detections. In Section 12.41 we discuss the sensitivity of the different systems for 
detecting the correlation and compare the detected systems with their sensitivity. 
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2.1. Transit timings 


We used the transit timing catalog obtained by 


Holczer et al. 


f)2015. in preparationi) . 


publicly available at: ftp://wise-ftp.tau.ac.il/pub/tauttv/TTV/ver_112, The catalog 
includes all the KOIs listed in the NASA Exoplanet Archivqj, except KOIs 


• labeled as false positives, 

• with orbital periods longer than 300 days, 

• with transit depth larger than 10 %, 

• with phase folded transit signal-to-noise ratio (SNR) below 7.1, 

• suspected to be binaries because of a difference in transit depth between odd and even 
transits, 

• with folded light curve that did not display a transit with high enough signihcance; 
the p-value of the J^-test between the transit model and the no-variability model (flat 
line) did not yield a level of signihcance higher than 10“^. 


After removing all these KOIs, we are left with 2,600 systems. 


The full list includes a total of 295,373 individual transits. Of those, 71,240 were 
rejected based on (1) an J^-test comparing the individual transit model with a hat light 
curve (no variability), (2) a strong deviation of the transit depth, duration, or timing 
residual (after subtracting a linear ephemeris) from the sample mean, or (3) the transit 
is too close to, or overlapping with , a transit of another planet in the same system (see 


Holczer et al. 


(12015. in preparation!) for details). This left 224,133 entries. 


^ http: / / exoplanet archive. ipac .caltech.edu / 
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2.2. The local photometric slope 


We derived the local photometric slope at each transit by htting a polynomial to the 
stellar light cnrve, extending fonr transit dnrations centered on the transit, while ignoring 
the in-transit data. For each transit event we normalized the ont-of-transit data to the 
median brightness of the relevant Kepler qnarter. 


The fitting was done in two steps. First, we fitted independently six polynomials, of 
degrees of one throngh six, using the MATLAB/REGRESSION function that exercises the 
linear least squares approach. The initial £t was followed by a hnal fit, after outliers beyond 
four times the scattein were ignored. In the second step we chose the best polynomial fit for 
each transit event by performing J^-tests between all pairs of polynomial fits and calculating 
their p-values. We chose the best fit to be the one with the highest degree for which the 
p-values of all the J^-test pairs with polynomial fits of lower degrees were lower than 10“^. 
To obtain the local brightness slope during the transit we calculated the derivative of the 
best-fit polynomial at mid-transit time. The derivative uncertainty was propagated from 
the errors of the polynomial coefficients. 


Next, we examined the sample of local slopes for each KOI and removed outlier slopes. 
A slope was identified as an outlier if it deviated from the sample mean by more than five 
times the sample scatter plus three times the uncertainties median. Among the 224,133 
measured local slopes we identified 265 (= 0.12 %) as outliers, leaving 223,868 to be used 
in the following analysis. 

Table [1] lists all transit events we analyzed, including the fitted local slope, s, the local 
slope uncertainty, ag, the degree of the best-£t polynomial, and a flag marking outlier local 


^Throughout this work we define the scatter as 1.4826 times the median absolute deviation 
(MAD), which equals the standard deviation for a Gaussian distribution. 





slopes. 


Table 1: Light curve local slope for KOI transits 


KOI 

n“ 

t ^ 

[d] 

[ppm/d] 

[ppm/d] 

deg® 

flag'" 

1.01 

0 

55.7633 

-39 

66 

1 

0 

1.01 

1 

58.2340 

-47 

81 

1 

0 

1.01 

2 

60.7046 

132 

68 

1 

0 

1.01 

4 

65.6458 

19 

70 

1 

0 

1.01 

5 

68.1164 

-40 

50 

1 

0 

1.01 

6 

70.5870 

-59 

83 

1 

0 

1.01 

7 

73.0576 

-6 

72 

1 

0 

1.01 

8 

75.5282 

14 

75 

1 

0 

1.01 

9 

77.9989 

17 

79 

1 

0 

1.01 

10 

80.4695 

-19 

58 

1 

0 


(This table is available in its entirety in a machine-readable form in ftp://wise-ftp.tau.ac.il/pub/tauttv/Slope 
A portion of the table is shown here for guidance regarding its form and content.) 


“Transit number. 

^ Expected transit t ime of the linear ephemeris in BJD - 
(12015. in preparationl) . 

^Derived local slope. 

'^Derived local slope uncertainty. 

“Polynomial degree of the best fit chosen. 

•^Transit flag: 1 = transit identified as an outlier. 0 = not an outlier. 


2454900 


following 


Holczer et al. 


2.3. Searching for slope-TTV (anti)correlation 


We now turn to identify KOIs with a statistically signihcant correlation between the 
TTV and the local photometric slope. We chose to derive the correlation by htting a linear 
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function to the TTV vs. the local slope (using the MATLAB/REGRESSION function), and 
then identify the KOIs for which the slope of that linear function strongly deviates from 
zero. 


We referred to the fitted linear function’s slope as the global slope, QS, with an 
uncertainty agg. We used the absolute value of the ratio of the two, QSjogs, as a proxy 
to the correlation’s statistical signihcance. We performed this £t only for KOIs with 20 
or more measured transits, to avoid the effects of small number statistics. Out of 2,600 
KOIs, 1,858 had 20 or more measurements. To assign a false alarm probability (FAP) to 
each correlation, we performed 10^ bootstrap tests, randomly permnting the sets of local 
slopes. We defined the p-valne of the derived global slope as the fraction of permntations 
that yielded QSjogs larger in its absolnte valne than that of the real data set. 


Table [2] lists the statistical parameters for the KOIs in the sample. I t also inclndes 
the orbital period, tr ansit depth, and the TTV uncertainties, taken from iHolczer et ah 


(12015. in preparation 


, and the nnmber of transits for which both the local slope and the 


TTV are not rejected as ontliers. The fitted QS, ogs, and the corresponding p-value are 
listed for 1,858 KOIs which had 20 or more measnrements. 


Our search could be applied only to active stars with 
modu l ation. We expected those stars to be identihed by 


(12018 


arge enongh spot-indnced ste llar 


McQuillan. Mazeh &: Aigrain 


2014I) . who searched the Kepler light curves for stellar rotation. Therefore while our 


analysis (as described above) was applied to all 2,600 KOIs, we have applied a more detailed 
analysis (as described below) to the snbsample of 862 KOIs whose host star rotation period 


and p hotometric activity amplitnde were measured by 


McQuillan. Mazeh fc AigrainI (12013 


20141 1. For KOIs in that snbsample we have calculated the maximum expected TTV due 


^ftp://wise-ftp.tau.ac.il/pub/tauttv/TTV /ver_112 


















to spot-crossing, max{TTV'sc}, by Equation (13) of Paper I, using the published stellar 
rotation period and amplitude, all listed in Table [2J Of those 862 KOIs, only 726 had at 
least 20 measurements of the TTV and the local slope. 


2.4. Detection sensitivity of the slope-TTV (anti)correlation 

We now turn to discuss the sensitivity of each KOI’s data to the detection of the 
TTV-local-slope correlation, assuming the time variation is induced by spot-crossing 
events. The detection depends on the magnitude of the induced TTVs, in terms of the 
TTV uncertainties, and on the square root of the number of measurements (i.e., transits). 
Therefore we dehne a new quantity, S, to be: 

5^max{rrK.4^_ (1) 

CtTV 

where axrv is the median TTV uncertainty and N the number of transit timing 
measurements. We note that we can derive S only for the 862 systems for which the 
photometric modulation of the stellar rotation was obtained. For those systems, all four 
quantities in Equation ([1]) are listed in Table |2l 

In Figured] we plot S against the derived p-value of the derived slope in logarithmic 
scale. Only the 726 KOIs with measured stellar rotation and more than 20 transits 
are plotted, as only those allow to calculate S. The upper (lower) panel shows KOIs 
with positive (negative) correlations. The dashed line in both panels is at 5=50, chosen 
somewhat arbitrarily as a threshold to identify KOIs with a strong sensitivity to TTV 
induced by spot-crossing. 

In order to single out the systems with significant correlation detection we picked out 
the KOIs that have a p-value lower than 10“®. We found six such KOIs with negative QS, 
suggesting a prograde conhguration, while not a single system with positive QS. In fact, we 
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Fig. 1.— Sensitivity to the correlation of TTV indnced by spot-crossing with the brightness 
slope, S, vs. the derived statistical significance of the correlation for the sample of 726 KOIs 
with detected stellar rotation and at least 20 transits. The npper panel shows KOIs for 
which we obtained a positive slope and the lower panel shows KOIs with a negative slope. 
The red crosses are the KOIs with p-valnes lower than 10“® and the green dashed line marks 
an expected TTV sensitivity of iS = 50. The latter is a somewhat arbitrary valne, used as a 
threshold to mark KOIs with a strong TTV sensitivity. 

did not hnd any system with a positive QS and a p-value below 10“^. These six KOIs with 
signihcant correlation are marked in the hgure by red crosses, and are discussed in detail in 
Section 01 
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One would expect to find systems showing a strong correlation (or anti-correlation) 
due to spot-crossings in the upper left corner, with high sensitivity and low p-value. Indeed, 
we hnd the six KOIs with signihcant correlation in that corner in the lower panel, while the 
same corner in the upper panel is empty. For hve out of these systems all our 10^ random 
permutations gave a QS/ags smaller (in absolute value) than that of the real data, so 
somewhat arbitrarily we assigned them a p-value of 5-10“® in Figured] We also hnd three 
other KOIs with high sensitivity (above the dashed line in both panels of Figured]) but with 
no signihcant detection. The nine systems and their parameters are listed in Table [S] This 
includes the stellar ehective temperature Teg, taken from the NASA Exoplanet Archive on 
December 1st, 2014, and the impact parameter b of the planetary orbit, derived here as 
described in Section [S] 


3. The impact parameter and its impact on the correlation 


As discussed and demonstrated in the simulations of Paper I, the TTV local-slope 
correlation can be detected only if the spot crossing events occur in varying locations on 
the stellar disc. This can happen only if the chords of the planet and the spot lie across 
the stellar disc side by side. If the latitude of the spot, 9, is such that cos^ is substantially 
different from the impact parameter of the transiting planet, or if the stellar obliquity is 
substantially different from zero (or from 180°), there cannot be any observed correlation. 
It is therefore of interest to examine the impact parameters of the systems with high 
sensitivity and see if the information about the planet motion over the stellar disc can 
tell us something about the location of the stellar spots. This discussion does not include 
KOI-1546, which was recently identihed to be an eclipsing stellar binarjc. 


^The true nature of this object was noted on March 2015 on 


the Community Follow-up Observing Program (CFOP) website. 


see: 
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We decided not to use the NASA Exoplanet Archive catalog values for the impact 
parameters of the KOIs, as they are based on long cadence data only, whereas short cadence 
data are available for all our high sensitivity systems. Furthermore, the KOI sample on the 
NASA Exoplanet Archive displays an impact parameter distribution with a strong increase 
towards 6 = 0 (at the time of writing), instead of the expected uniform distribution. We 
do not know the cause for this artifact, but it suggests the use of those values may be 
unreliable for our tests. 


Therefore, we derived the impact parameters ourselves for the eight systems, using the 
short cadence simple-aperture photometry data. The hrst step was detrending the transit 
light curve by Ending the continuum around each transit, ignoring the photometry of the 
transit itself, up to 0.6 transit durations around the expected transit center. We then fitted 
a more extended region, of up to four durations around the expected transit center, and 
fitted six different polynomials of degrees one to six. An J^-test was performed between 
all of the different pairs of polynomial fits. The best fit was chosen as the one with the 
highest degree, for which the p-value of all the J^-tests with polynomial fits of lower degrees 
was lower than 10“^. The last stage of the detrending process was to add back the transit 
points and divide the flux at that region by the chosen polynomial fit. 


To obtain a model for the transit we folded the light curves of the different transits 


together using the 


oest-fitting period an d phase. We then modeled the transit using the 


OCCULTNL model flMandel fc Ago 


20021) . where the position of the planet was chosen 


with a rectilinear model of normalized impact parameter b, transit duration T^ur from 
mid-ingress to mid-egr ess, a linear limb darkening law with parameter C 2 (and other c* of 


Mandel fc Agoll (120021) fixed to zero), and planet-to-star radius ratio Rp/R-^. The model 


assumed no third light for the eight systems, consistent with the small contamination 


https://cfop.ipac.caltech.edu/edit_target.php?id=1546 
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(< 


10%) reported at the MAST. 


The ht was performed via the non-linear fitter IDL/mpfit flMarkwardti 


200911 ■ The 


scatter of out-of-transit points, a, was used to determine the uncertainty of the individual 
measurements, instead of relying on the reported photometric errors. Ten iterations of the 
fitter were run, followed by rejection of points lying more than lOcr away from the model; 
this procedure was repeated; and a hnal ten iterations of the htter resulted in the impact 
parameters reported in Table |3] and plotted in Figure |2j The reduced varied from 0.99 
to 1.78 with a mean of 1.23—these moderately exceed what white noise would have given, 
and they may be due to the scatter induced by the star spot crossing events. 

Figure |2] shows that the impact parameters of the five systems displaying TTV 
correlated with the local slope are in the range of 0.2 ^ b < 0.5. Interestingly, the Sun’s 
spots reach up to about 30° in latitude, and their size grows smaller as they descend in 
l ing the equato r when the Sun’s activity cycle is at minimum (e.g.. 


latitude unti 

reac 

Wilson et al. 

1996 


Li et ah 


2000l ). Therefore, to an observer viewing the Sun from the 


equator plane, transiting planets would cross sunspots only if they have b < 0.5, and would 
be difficult to detect at 5 ~ 0, just as we see in these systems. 

The only system with high sensitivity and b < 0.5 that does not show detected 


correlation is KOI-3.01 (= Kepler-3b = HAT-P- llb: 


Hirano et al. 


2011 


Sanchis-Oieda fc Winn 


2011 


Bakos et al. 


; see also 


2010 


Bekv et al. 


Winn et al 


2010b 


2014). This may be 


accounted for by the large stellar obliquity of KOI-3.01, as explained in more details in 
Section O A si milar scenario might acco unt for the lack of a correlation f or KOI-63.01 


= K epler-63b; 


Sanchis-Oieda et al. 


2oi3l . For KOI-254.01 (= Kepler-45b; 


Johnsonetah 


2 OI 2 I ) there is no evidence for correlation and therefore for spot crossings, which might 
be explained by the relatively large impact parameter, beyond 0.5, that causes the planet 
to move across the stellar disk at higher latitudes than the spots. These three systems 
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Fig. 2.— Sensitivity to the correlation of TTV indnced by spot-crossing with the brightness 
slope, S, vs. the derived impact parameter of the ir transiting planet. Impact parameter 


was derived from the short cadence data, using the 


Mandel fc Ago! (120021) model. For some 


systems the error bar of the impact parameter is smaller than the marker size. The red 
points are the KOIs with detected correlation, with p-value lower than 10“^, and the blue 
points are the systems with no signihcant detection of a correlation. 


(KOI-3.01, KOI-63.01, and KOI-254.01) are reviewed in more details in Section [5l 

Taken together, although extremely small, our sample suggests that the stars considered 
here have spots within 30° of their equator but not exactly along the equator, in analogy 
to sunspots. However, this does not exclude the existence of high latitude star spots (e.g.. 
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Li et al. 


20001) since our method is less sensitive to such systems. This is so because as 


the spot moves to higher latitude both the stellar photometric modulation amplitude and 
transit timing precision decrease. 


In the next two sections we present the detailed analysis of the KOIs with high 
sensitivity. Section 0] discusses the systems with signihcant correlation, and Section O 
discusses the systems with non-signihcant correlation. 


4. The six KOIs with significant correlation 


In this section we present the six systems that show a highly signihcant anticorrelation 
between their TTV and the corresponding local slope, with QSjags p-value lower than 
10“^. Their analysis, based on the long cadence, are presented in Figures [3HTH Figures [3l 
El [3 [H [m and [l3] show the light curve analysis, and Figures HI El El EQI [l2l and EH show the 
TTV analysis. 

In Figures [3l El 171 [9l ET] and EH the left panel presents the phase folded light curve 
around the transit, overplotted by the htted model (red dashed line), with the residuals at 
the bottom. In all six systems, the larger scatter in the residuals during transit compared 
to out of transit is probably due to spot-crossing events. 


The top right panel presents a segment of the light curve as a function of time, 
where the photometric modulation due to stellar rotation and activity is clearly noticeable. 
We can see that the modulation is not strictly periodic, as the spots evolve in size and 


loc ation. The spacing between t 


by 


l e vert i cal re d dashed lines is the rotation period identihed 


McQuillan. Mazeh fc Aigrain 


(120131, 


20141) . They used the Auto-Correlation Function 


(ACF), presented in the bottom right panel, where the red dashed line marks the rotation 
period, which, as expected, coincides with the shortest time lag ACF peak. 
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In Figures HI |6l [H [TOl [12] and [HI the top left panel shows the TTV vs. the local 
photometric slope with the linear ht overplotted as a solid red line. The clear anti-correlation 
between the TTV and the local slope indicates prograde motion. The top right panel shows 
the measured TTV as a function of time (top) and the corrected TTV (bottom) after 
subtracting the linear £t of the TTV vs. local slope. The corrected TTV does not show 
much difference relative to the measured TTV (except maybe in KOI-1546.01). 


The bottom left panel shows the TTV frequency power spectrum, with a green dashed 
line marking the stellar rotation frequency. When relevant, another green dotted-dashed 
line or two mark the harmonics of this frequency. The red dotted-dashed line marks 


the sampling stroboscopic 


frequency—^^the T 


Kepler long cadence data fjSzabo et ah 


2013 


'V signal induced 


Mazeh et ah 


by the sampling of the 


201311 . The power spectrum is 


plotted up to the Nyquist frequency, which in this case is half the orbital frequency, at 
which the TTV signal is sampled. The bottom right panel shows the TTV ACF with the 
stellar rotation period marked by a green dashed line. 


In two cases, KOI-203 and KOI-883, we see a clear TTV periodicity with the rotational 
period. This is another indication that the TTV signal is related to the stellar rotation. 

In the case of KOI-217.01 and KOI-1074 1 Figures [61 and [TTH the strongest power spectrum 
(PS) peak is at twice the stellar rotation frequency. This is consistent with a scenario in 
which the star has two significant spots (or groups of spots), each at opposite hemispheres, 
as can be seen also in the stellar light curves (Figures [5] and [13] top right panel). 

Note that for KOI-895, KOI-1074, and KOI-1546.01, the rotational frequency /rot is 
larger than the Nyquist frequency of the sampling, which is half the orbital frequency, 
/xyq = 0.5/orb (scc the Orbital and rotation periods listed in Table [3]). Therefore the 
rotational frequency is aliased to a smaller frequency. 

Out of the six KOIs, two are confirmed planet hosts. KOI-203.01 (Figure [3]) was 
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confirmed by 


Desert et al. 


(12011 


, Kepler-17b), who used Kepler short cadence data and 


found it to be prograde with a spin-orbit angle smalle r than 15°, which a grees well with our 


hndings. The second conhrmed planet is KOI-217.01 flHowell et ah 


20101. Kepler-71b). 


5. KOIs with no significant correlation 


In this section we present the three KOIs with high sensitivity but with no detected 


correlation between the TTV and the corr esponding local-slope. All t 


iree. KOI-3.01 (= 


Kepler-db = HAT-P-llb: Bakos et al. 

>nini. KOi-63.n 

1(=] 

2013[). and KOI-254.01 (= Kepler-45b: 

Johnson et al. 

2012) 


Sanchis-Oieda et ah 


2 OI 2 I) are confirmed planets. They 


are presented in Figures [THl through 1201 in the same format as the KOIs showing significant 
correlation (see Section H] and Figures 131 through [TT|) . It can be visually seen that the TTV 
vs. local slope does not show a signihcant correlation for all three systems (see Figures fT6l 

[T51 and [20]). 


In order to allow a more careful examination of the light curves of these three systems, 
and specihcally to try and understand why they do not show a s ignific ant correlation, we 


use here the Kepler short cadence data and the 
generating Figures [15] through |20l 


Mandel fc Agoll ( 20021 ) transit model for 


For KOI-3.01 (= Kepler-3b = HAT-P-llb; see Figures [Thl and [TOl) the residuals of the 
phase folded and binned light curve are not consistent with white noise but instead show 
clear systematic features (Figure [15] left panel). This is in contrast to the other systems 
(see Figures [3l [5] [7] [H [TTl and[T3|) where the residuals of the in-transit light curve show an 
increased noise level but no systematic features. This distribution of the residuals suggests 
that the same features occur in many transit events at the same phases, so they do not 
average out in the folded light curve. Assuming these residuals are caused by spot-crossing 
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events, this might suggest that the spot crossings occur primarily at specific phases during 
the transit. 


Fort unately, the Keyler light curve of KOI-3.01 was 


authors fjSanchis-Oieda fc Winn 


observations of the RM effect (jWinn et ah 


2011 


Demine et ahl 


2ninbl: 


2011 


carefully studied by several 


Beky et ah 


Hirano et ah 


2014) . in addition to 


201ll) . These studies found 


that the system has a sky-projected spin-orbit angle close to 90°, and that the host star 
has two active latitudes. Hence, spot-crossing events occur when the planet moves across 
the active stellar latitudes, explaining why they are seen only at specific phases during the 
transit. Since those phases sample only part of the transit and not the entire transit it 
reduces the TTV range induced by the spot-crossing events. This is also consistent with a 
weak TTV periodicity identified at the rotation period (see Figure [T6|) . where no significant 
peak is seen at the rotation frequency (bottom left panel) and only a small peak, close to 
the noise level at the ACF (bottom right panel). Moreover, in this configuration, when a 
spot is crossed by the planet the star is at a similar phase in its rotation, resulting in a 
small range of local slopes. The combined result is that although this system posses the 
sensitivity for TTV-local-slope correlation, its characteristics suppress that correlation. 

KOI-63.01 (= Kepler-63b; see Figures fTTI and flSD also s hows systematic feat u res in 


the phase folded light curve residuals (Figure ITTI left panel). 


Sanchis-Ojeda et al. 


(120131) 


have already studied the Kepler-63b Kepler light curve and observed the RM effect. They 
determined the system has a sky-projected spin-orbit misalignment of - 1101^4 
a stellar spin axis inclination angle of 138±7 deg. Therefore, in this configuration the 
spot-crossing events are confined to the same parts of the transit. This explains the features 
in the phased light curve residuals and the non-detection of a spots-induced TTV signal, 
where the latter leads to no TTV-local-slope correlation. Yet, the TTV period analysis 
(Figure fTSjl does show a significant TTV periodicity at the rotation frequency, which in this 
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case is aliased to /rot,alias = 2/orb - /rot, as 3/Nyq < /rot < 4/Nyq. This could be explained by 
assuming there are one or two spot latitudes and therefore one or two transit phases that 
allow spot crossing, not enough to produce a detectable correlation, but enough to induce a 
TTV modulation with the rotation period. 


For KOI-254.01 (= Kepler-45b, 


Johnson et ah 


2012 


; see Figures IT^ and [20]) . the 


spin-orbit angle is currently not known. We hypothesize that the planet does not cross 
spots during the traii sit and that the system is spin-orbit aligned, similar to, e.g., Kepler-77 


(iGandolh et ah 


201311 . This is consistent with the noise level of the in-transit light curve 


being similar to the noise level of the out-of-transit phases, which is different from the six 
systems with signihcant correlation that show an increased in-transit residuals noise level 
(see Section 0] and Figures |3l |5l [71 |9l [TTl and [13]). For this syste m we have determin ed here 


a transit impact parameter of 0.5732 ± 0.0027, consistent with 


■Tohnson et ah 


2012 


value 


of 0.6 ± 0.2, who used less Kepler data than used here. This impact parameter is larger 
than the ones of the systems with signihcant correlation (see Figure [2]), suggesting that 
KOI-254.01 impact parameter is larger than the impact parameter corresponding to the 
typical active stellar latitudes. For this object the TTV period analysis (Figure [20]) does 
not show a periodicity at the rotation period, as expected. 


Another possible scenario that can weaken the correlation between the TTV and 
the local slope and was not mentioned above is a dynamical TT V signal, following th e 


gravitational interaction with anot 


Ford et ah 

2012; 

Steffen et ah 

2012 


ler body in the systems (e.g.. 


Fabrvckv et ah 


2012 


20121) . However, there is no evidence for such a signal in any 


of the three systems discussed above. 


To summarize this section, for all three systems we are able to provide plausible 
scenarios that explain why despite having the sensitivity to a TTV signal induced by 
spot-crossing events, these systems do not show a detected correlation. 
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6. Discussion 


We presented here the application of a simple method to distinguish between prograde 
and retrograde transiting star-planet systems, where the brightness of the host star is 
modulated by stellar spots, and the transits include spot-crossing events. The method is 
based on the assumption that even when spot-crossing events are not resolved, they induce a 
shift in the derived mid-transit timing, when fitting a model that ignores the spot-crossings. 


We have applied o ur method to a samp 


by 


Holczer et al 


e of 2,600 KOIs whose transit timings were measured 


(12015. in preparation!) using the long cadence data, and are publicly 


available (ftp://wise-ftp.tan.ac.il/pub/tauttv/T' 


systems with published stellar rotation periods flMcQuillan. Mazeh fc Aigrain 


'V/ver_112). We have concentrated on 86 2 


2013 


20i4 . 


out of which only 726 had at least 20 measurements of the TTV and local slope. 


Using the formalism of Paper I, we identified nine systems, listed in Table 3, with high 
enough stellar modulations and large enough planets to allow the detection of the TTV 
local-slope correlation. The nine systems have relatively large photometric modulation 
amplitudes due to rotation (> 5000 ppm), all of them well above the median of the whole 
analyzed sample (~3000 ppm). The precision of the transit timing for these nine systems 
is high, with a median of 0.45 min, while the sample median is ~13 min. This indicates 
that our method requires objects showing relatively strong rotation modulation and high 
SNR transits, the latter leading to highly precise mid-transit timings. The nine systems we 
have identified have orbital periods in the range of 0.9-9.4 d and transit depth of 0.4-4.0%, 
consistent with the typical parameters range of hot Jupiters. The small size of the sample 
is not surprising, given the fact that we are looking for a minor effect that is hiding in the 
noise in most cases. One of the nine cases is KOI-1546, which recently has been found to be 
a stellar binary, so we are left with eight bona fide planet-candidate systems (five of which 
already confirmed as planets). 
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It turned out that out of these eight systems five show a clear negative correlation, 
indicating prograde motion, and none a positive correlation, which could hav e indicated a 
retrogr ade motion. This is consistent with the hnding of the seminal work of 


Winn et ah 


(j2010aj), who showed that most cool planet-host stars have close to zero obliquity, as four 


out of our hve systems have lower surface temperature. However, this could also result from 
an intrinsic bias in our approach, since it requires active host stars, which are typically 
convective and hence relatively cold. 


As shown in Figure |21 the hve systems with detected correlation (not including 
KOI-1546, which is a binary) all have an impact parameter in the range of 0.2 < 6 < 0.5. 
Assuming those systems have their spots moving on similar latitudes as the planetary chord 
as seen from Earth, this result shows a similarity to the Sun, where the spots reach up to 
about 30° latitude and decrease in the size as they move closer to the Sun’s equator. 


To summarize, we have found evidence for prograde motion of transiting planets 
but no retrograde motion. Furthermore, we found some hints that stellar spots might 
preferentially be at low, Sun-like, latitudes on the stellar disc (< 30°) but not exactly along 
the stellar equator, although our search suffers from an observational bias against detection 
of high latitude spots, and the detected sample is extremely small. The original goal of this 
project was to learn about the planetary motion, but we were surprised to hud a possible 
hint about a feature of spot motion. In principle this is possible, as the planet with its 
small radius is sampling the stellar surface with high resolution that is not available in 
present-day observational techniques. The high precision, evenly spaced, long-span light 
curves of Kepler make this possible. 


The approach described here is not limited to K epler star-p 


principle, be applied to Kepler stellar binary systems fjPrsa et ah 


anet systems and can, in 


2011 


Slawson et ah 


201ll), 


an analysis that we are preparing in a forthcoming paper. In fact, a negative correlation 
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between TTV and local slope was identified already for the eclipsing ste llar binary system 


within the circumbinary planetary system Kepler-47 fjOrosz et al 


2 OI 2 I), where a rnore 


detailed analysis of the spot-crossing events indicates a prograde motion flOrosz et al 


201211 ■ 


In the future the method used here can be applied to a large sample of systems 


monitored by current and 


inclu ding K2 flHowell et al 


hture space-ba sed surveys c 


2ni4ll . TESS flRicker et ah 


ehyering high-qualit y photometry , 


2014 1. and PLATO f Rauer et al 


2014fl . A larger sample can reveal the typical characteristics of prograde and retrograde 
systems, which in turn will allow constraining formation and orbital evolution processes. 
In addition, a larger sample of systems where spot-crossing events are identified will also 
better constrain spots’ behavior. 
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Table 2: Statistical parameters of the TTVs and the stellar rotation of the KOIs 


KOI 

Period^ 

Depth'^ 

Period^ 

Amp.'^ 

flag® 

QS^ 

crgs^ 

p-iGS/ags)^ 

Expected’ 

(JTTV^ 

Nk 

5' 


Orb. 


Rot. 

Rot. 





TTV 





[d] 

[ppm] 

[d] 

[ppm] 


[min d] 

[min d] 

[log] 

[min] 

[min] 



I.Ol 

2.47 

14210 



0 

-46 

61 

0.0 


0.08 

428 


2.01 

2.20 

6694 



0 

-32 

76 

-0.1 


0.23 

599 


3.01 

4.89 

4361 

29.47 

5875 

1 

26 

19 

-0.7 

2.97 

0.24 

214 

178.8 

5.01 

4.78 

980 



0 

-300 

1000 

0.0 


1.65 

278 


7.01 

3.21 

736 



0 

1100 

1900 

-0.3 


3.53 

325 


10.01 

3.52 

9370 



0 

130 

no 

-0.6 


0.60 

380 


12.01 

17.86 

9228 

1.25 

390 

1 

70 

150 

-0.2 

0.45 

0.37 

71 

10.3 

13.01 

1.76 

4602 



0 

196 

62 

-4.7 


0.13 

735 


17.01 

3.23 

10811 



0 

50 

120 

-0.7 


0.33 

334 


18.01 

3.55 

7454 



0 

1060 

190 

-7.0 


0.57 

376 



(This table is available in its entirety in a machine-readable form in ftp://wise-ftp.tau.ac.il/pub/tauttv/Slope 
A portion is shown here for guidance regarding its form and content.) 


“Orbital period. 
^Transit depth. 
^Rotation period. 



and no period was identified, 1: A rotation period was detected, 2: KOI was not analyzed. 
■^Global slope. The slope of the linear fit to the TTV as a function of the local slope. 

® Global slope uncertainty. 

^Log p-value of QS/ags- 
“The maximum expected TTV. 

•^The TTV uncertainties median. 

^Number of measurements. 

^TTV Sensitivity. 
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Table 3: Statistical parameters of the KOIs with high TTV sensitivity 


KOI 

Period‘s 

Depth^ 

Period'^ 

Amp."^ 



p-(gS/ags)<^ 

Expected** 

(^TTV^ 

Nl 


Teff* 



Orb. 


Rot. 


Rot. 



TTV 







[d] 

[ppm] 

[d] 

[ppm] 

[min d] 

[min d] 

[log] 

[min] 

[min] 



[K] 


3.01 

4.89 

4361 

29.47 

5875 

26 

19 

-0.7 

2.97 

0.24 

214 

178.8 

4777 

0.28 

63.01 

9.43 

4037 

5.41 

8910 

2 

10 

-0.1 

4.27 

0.44 

139 

114.9 

5650 

0.73 

203.01 

1.49 

21575 

12.16 

12865 

-33 

3 

r- 

1 

V 

3.20 

0.29 

710 

292.1 

5624 

0.24 

217.01 

3.91 

22620 

19.77 

4970 

-140 

20 

< -7 

1.46 

0.45 

331 

58.4 

5543 

0.26 

254.01 

2.46 

40367 

15.81 

10930 

-9 

6 

-0.8 

1.57 

0.48 

540 

75.4 

3820 

0.57 

883.01 

2.69 

39186 

9.02 

11340 

-38 

4 

< -7 

1.94 

0.40 

495 

107.1 

4809 

0.51 

895.01 

4.41 

13822 

5.07 

10575 

-39 

5 

A 

1 

5.33 

0.95 

302 

97.0 

5600 

0.37 

1074.01 

3.77 

14018 

4.04 

8360 

-32 

6 

-5.6 

4.12 

1.04 

296 

67.9 

6302 

0.36 

1546.01 

0.92 

15752 

0.91 

7600 

-27 

1 

1 

V 

1.90 

0.62 

1430 

115.6 

5713 

0.67 


“Orbital period. 
^Transit depth. 
'^Rotation period. 


‘^Rotation semi amplitude (see 


McQuillan. Mazeh fc Aigrain 


2013 for details). 


“Global slope. The slope of the linear fit to the TTV as a function of the local slope, 
t Global slope uncertainty. 

®Log p-value of GS/ags- 
^The maximum expected TTV. 

*The TTV uncertainties median. 


•^Number of measurements. 

'^TTV Sensitivity. 

^Stellar effective temperature taken from the NASA Exoplanet Archive. Typical uncertainty is 100-200 K. 
Impact parameter derived here using Kepler short cadence data (see Section |3] below). Typical uncertainty 
is 0.01-0.02. For KOT1546.01 no short cadence data is available so we used the value from the NASA 
Exoplanet Archive. 
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Fig. 3.— KOI-203.01 (= Kepler-17b): Transit light curve and stellar modulation and its 
autocorrelation. Left: Model fit to the detrended phase folded transit light curve (top) 
and its residuals (bottom, in ppm). The htted model is overplotted in a red dashed line. 
Upper right: A segment of the light curve. The spacing between the vertical red dashed 
lines equals the host star’s rotation period. One can follow the stellar rotation modulation 
with the detected period. The modulation changes its shape and amplitude, but retains its 
period. Bottom right: The Auto-Correlation Function (ACF), with a vertical red dashed 
line where the rotation period is found. One can clearly see the ACF peaks at lags equal to 
integer multiples of the rotation period. 
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Fig. 4.— Analysis of the TTV of KOI-203.01 (= Kepler-17b). Upper left: The TTV, in 
minntes, as a fnnction of the corresponding local slope, in ppm per thousand days. The 
overplotted solid line (red) is a linear £t whose slope is defined as the global slope. Upper 
right: The derived TTVs as function of time. The upper panel shows the derived TTV, 
while the lower panel displays the corrected one, after subtracting the linear £t plotted in 
upper left panel. Bottom left: The TTV power spectrum. One can clearly see strong peaks 
corresponding to the stellar rotation period, marked by a green dashed line, and its hrst 
two harmonics, marked by green dotted-dashed lines. The red dotted-dashed line marks the 
sampling stroboscopic freq uency which is the TTV signal induc ed by the sampling of the 


Kepler long cadence data fISzabo et al 


2013 


Mazeh et al 


20131) . The power spectrum is 


plotted up to the Nyquist frequency which in this case is half the orbital frequency. Bottom 
right: The ACF of the derived TTV. The green dashed line presents the lag corresponding 
to the stellar rotation. One can notice a small peak at this lag and another one at a lag 
twice the rotation oeriod. 
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Fig. 5.— KOI-217.01 (= Kepler-71b) transit light curve and stellar modulation and its 
autocorrelation. See Figure [3] for details. 
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Fig. 6.— Analysis of the TTV of KOI-217.01 (= Kepler-71b). See Figure IHfor details. The 
green dashed line in the bottom left panel mark the stellar rotation frequency, and its hrst 
harmonic is marked by a green dotted-dashed line. 
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Fig. 7.— KOI-883.01 transit light curve and stellar modulation and its autocorrelation. See 
Figure [3] for details. 
















































Power Spectrum TTV [min] 



-20 -10 0 10 20 
Local Slope [ppm/1 OOOd] 






! 





1 






















0.05 0.1 0.15 

Frequency [1/d] 


2 

0 

„-2 

c 

^-4 

> 

1= 2 
0 

-2 

-4 


200 400 600 800 1000 1200 1400 

Time [d] 



Lag [d] 



Fig. 8.— Analysis of the TTV of KOI-883.01. See Figure H] for details. The green dashed 
line in the bottom left panel marks the stellar rotation frequency. 
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Fig. 9.— KOI-895.01 transit light curve and stellar modulation and its autocorrelation. See 
Figure [3] for details. 
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Fig. 10.— Analysis of the TTV of KOI-895.01. See Figure H] for details. The green dashed 
line in the bottom left panel mark the stellar rotation frequency, and its hrst harmonic is 
marked by a green dotted-dashed line. 
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Fig. 11.— KOI-1546.01 transit light curve and stellar modulation and its autocorrelation. 
See Figure [3] for details. 
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Fig. 12.— Analysis of the TTV of KOI-1546.01. See Figure H] for details. The green dashed 
line in the bottom left panel marks the stellar rotation frequency. 

































































Res. [ppm] Rel. Flux 


41 




Time [d] 



Lag [d] 


Fig. 13.— KOI-1074.01 transit light curve and stellar modulation and its autocorrelation. 
See Figure 13] for details. 
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Fig. 14.— Analysis of the TTV of KOI-1074.01. See Figure 0] for details. The green dashed 
line in the bottom left panel mark the stellar rotation frequency, and its hrst harmonic is 
marked by a green dotted-dashed line. 
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Fig. 15.— KOI-3.01 transit light curve and stellar modulation and its autocorrelation. See 
Figure [3] for details. 



































Power Spectrum i i i TTV [min] 


44 



Frequency [1/d] 


Lag [d] 


Fig. 16.— Analysis of the TTV of KOI-3.01. See FigurelUfor details. The green dashed line 
in the bottom left panel mark the stellar rotation frequency. 
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Fig. 17.— KOI-63.01 transit light curve and stellar modulation and its autocorrelation. See 
Figure [3] for details. 
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Fig. 18.— Analysis of the TTV of KOI-63.01. See Figure H] for details. The green dashed 


line in the bottom left panel mark the stellar rotation frequency. 
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Fig. 19.— KOI-254.01 transit light curve and stellar modulation and its autocorrelation. 
See Figure [3] for details. 
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Fig. 20.— Analysis of the TTV of KOI-254.01. See Figure H] for details. The green dashed 
line in the bottom left panel mark the stellar rotation frequency. 



































































































































































































































































